Stabilization of the Electroweak Z String in the Early Universe 
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The standard electroweak theory admits a string solution, the Z string, in which only the electrically 
neutral Higgs fields are excited. This solution is unstable at zero temperature: Z strings decay by 
exciting charged Higgs modes. In the early Universe, however, there was a long period during which 
the Higgs particles were out of equilibrium but the photon field was in thermal equilibrium. We 
show that in this phase Z strings are stabilized by interactions of the charged Higgs modes with 
the photons. In a first temperature range immediately below the electroweak symmetry breaking 
scale, the stabilized embedded defects are symmetric in internal space (the charged scalar fields are 
not excited). There is a second critical temperature below which the stabilized embedded strings 
undergo a core phase transition and the charged scalar fields take on a nonvanishing value in the core 
of the strings. We show that stabilized embedded defects with an asymmetric core persist to very 
low temperatures. The stabilization mechanism discussed in this paper is a prototypical example of 
a process which will apply to a wider class of embedded defects in gauge theories. 



PACS numbers: 98.80.Cq, 11.27.+d 
1. Introduction 

In a previous paper [1] we suggested that a wide class 
of embedded defects may be stabilized by plasma effects. 
We studied a toy model consisting of four real scalar fields 
(f)i i — 1, ..,4 with the "generalized Mexican hat" poten- 
tial 



(1) 



two of which {(pi and (j>2 to be specific) being electrically 
charged, the other two neutral. At zero temperature, the 
vacuum manifold in this model is and hence does not 
admit any stable topological defects. However, there are 
embedded defects [2,3], configurations which span only 
a subset of the vacuum manifold. One example is the 
global cosmic string solution of the subspace of the the- 
ory with 01 = <?!'2 = 0. At zero temperature this string 
configuration is unstable and can decay through the ex- 
citation of the charged fields. In more graphic terms, 
the field configuration slips off the top of the potential in 
the charged field directions and approaches the vacuum 
manifold everywhere in space. However, in the presence 
of a bath of photons, interactions of the photons with 
the charged scalar fields will lead to an effective poten- 
tial which is lifted in the directions of the charged fields, 
thus generating a potential barrier which can stabilize 
the embedded string. This toy model is a theory with 
a global symmetry, and is realized in the Sigma model 
description of the low energy limit of Quantum Chro- 
modyanamics in the limit of vanishing pion mass (see 
[4] which includes a detailed discussion of effects which 
arise when the explicit symmetry breaking caused by the 
nonvanishing pion mass is taken into account). 

The standard electroweak theory is a good example 
of a theory which contains embedded strings in which 



only the neutral scalar field components of the multi- 
dimensional order parameter are excited. Among all pos- 
sible embedded strings, (see e.g. [5] for a classification 
of embedded strings in the electroweak theory) the elec- 
troweak Z string [6,7] is the string configuration consist- 
ing of excitations of the neutral fields only. 

In the paper we demonstrate that Z strings are stabi- 
lized during the temperature interval between the elec- 
troweak phase transition temperature and the recombi- 
nation. At a certain temperature during this period, a 
core phase transition takes place [8] below which the em- 
bedded strings will become superconducting as shown in 
[4] , thus enhancing their stability and leading to the for- 
mation of vortons [9]. The enhanced stability of elec- 
troweak strings due to current-carrying zero modes is 
similar to the stabilization by neutrino zero modes con- 
sidered in [10]. 

The mechanism discussed in this paper generalizes to 
a wide class of embedded defects. Since defects can play 
an important role in cosmology - even if they are not 
stable at all times - the mechanism discussed here may 
have important consequences for many aspects of cosmol- 
ogy. One possibility is to use embedded global anoma- 
lous strings such as the pion string of low energy QCD 
[II] to generate primordial magnetic fields [12]. Crucial 
aspects in this application are, besides the stabilization 
of the embedded strings, the fact that, via the anomaly, 
charged zero modes on the string generate coherent mag- 
netic fields circling the string [13], and the fact that the 
length scale of the string network increases in comoving 
coordinates (see e.g. [14]), which provides the mechanism 
for generating the required large coherence length of the 
primordial magnetic fields. QCD at large baryon den- 
sity also leads to the stabilization of a type of embedded 
strings, K-strings [15]. 

Stabilized embedded defects (in particular embedded 
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walls) could also play a role in defect-mediated baryo- 
genesis (see e.g. [16] and [17,18] for the main ideas of 
defect mediated GUT and electroweak baryogenesis, re- 
spectively), and they may be useful in implementing the 
scenario of Dvali et al. [19] for alleviating the monopole 
problem via defect interactions (see [20,21] for studies of 
the basic interaction mechanism). 

As with any class of topological defects, there will be 
severe cosmological constraints on models which admit 
them, resulting from the fact that the evolution of the 
defects in the early Universe may lead to predictions 
which arc in conflict with observations. The vorton abun- 
dance problem (sec e.g. [22]) leads to severe constraints, 
as does the constraint (specific to decaying defects) that 
the decay not lead to spectral distortions in the cosmic 
microwave background [23]. 

A new feature which arises in the present discussion of 
the electroweak theory, compared to the previous studies 
of toy models with a global symmetry, is the fact that 
the underlying symmetry is a local symmetry. 

2. Stabilization Mechanism 

Starting point is the Lagrangian for the standard elec- 
troweak theory: 



where the subscript NO stands for the Nielscn-Olesen 
U(l) cosmic string configuration [24], which in cylindrical 
coordinates r and 9 (the coordinates in the plane perpen- 
dicular to the string) and in temporal gauge A^^q = 
has the form 



$ivo(r,0) =r7./(^)e'^ 

A^,No{r,e) = - —(5^0; 
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where the indices fi and v are Lorentz indices, the index 
a is an SU(2) index, g and g' are the SU(2) and U(l) 
coupling constant, respectively, and W and Y are the 
SU(2) and U(l) field strength tensors, respectively. 

The field $ is a complex Higgs doublet, the upper com- 
ponent having positive charge, the lower component 
$° being neutral (if the expectation value of $ is chosen 
such that only the lower component is non- vanishing) . 
The analogy with (1) is clear: the lower SU(2) doublet 
component gives the two neutral real scalar fields 03 and 
04, the upper charged component yields 0i and 02- 

We are interested in describing the physics below the 
phase transition temperature Tc- In this case, the only 

gauge field which is excited is the photon field A^. In figuration reduces to the Z string), Z) is the Nielsen- 



where /(r) and v{r) arc the profile functions of the string. 
Note that since for the electroweak Z-string the vacuum 
expectation value of the upper component of $ always 
vanishes, the upper component of $ will be associated 
with electrically charged degrees of freedom everywhere 
in space. 

At this point the reader may object and recall that in 
the electroweak theory there is only one physical Higgs 
degree of freedom. At any point in space one can choose 
a gauge in which only one of the neutral Higgs fields 
is non-vanishing. The three scalar field degrees of free- 
dom corresponding to rotations in the vacuum manifold 
are eaten by the gauge fields acquiring masses, and it is 
only the massive Higgs degree of freedom which remains. 
Howevcir. the existence of topological defects in gauge 
theories is precisely a consequence of the fact that there 
are restrictions on the ability to impose the abovemen- 
tioned gauge uniformly over space. By causality, there 
will be regions in space in which the Higgs field is not 
in its vacuum manifold. This corresponds to localized 
energy configurations which cannot be gauged away. To 
describe the physics of such defects, it is advantageous to 
use a gauge in which the Higgs fields are not fixed, as we 
do below. 

In the absence of any excitations in the bulk, the elec- 
troweak Z string is unstable [25]. The unstable mode 
corresponds to an excitation of the charged Higgs dou- 
blet: 



$(C, x) ^ cos(C)$"(cos(0:c) sin(0$^ 



Z^(e,f)=cos(0^r(cos(e)x) 



(6) 



where ^ is the deformation parameter (for ^ = the con- 

(0) 



terms of the SU(2) gauge field W"" and the U(l) gauge 
field Y (we are suppressing the Lorentz index) , the Z and 
A fields are given by 



Z 



A = sm{e^)W^ + COS 



'■{0^)Y, 



(3) 



Olesen gauge field configuration, and <I>° is a complex 
Higgs doublet whose neutral component takes on the 
Nielsen-Olesen string configiuation, and whose charged 
component vanishes. The; field stand for the; complex 
Higgs doublet in which only the charged complex scalar 
is excited, and its transpose of is given by 



where 9^ is the weak mixing angle. 

The electroweak Z string is obtained by setting 



($+)r = ^(1^0) , 



(7) 



= W^ = A 
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Since by escaping into the charged scalar field direc- 
tions, the string configuration can decrease its poten- 
tial energy in the string core region, the energy per imit 
length E[$^) of the configuration (6) is smaller than the 
energy per unit length Eq of the Z string. Since the loss 
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in energy density is proportional to ^^y(0), and since the 
loss in energy density is confined to the core region of the 
string (core radius Vc), we obtain the following estimate 
for the maximal energy loss: 

E{0 >Eo- nXev'rl , (8) 

where k is a constant of order 1. A detailed analysis [25] 
shows that for realistic values of the weak mixing angle, 
there is indeed an energy loss. 

Another way to view the instability of the electroweak 
string at zero temperature [26] is as an instability to the 
formation of a W-condensate in the core of the string, a 
process which was shown to lower the energy. In [27] it 
was shown that this instability can be gauge transformed 
to an instability lowering the winding number N of the 
pure electroweak string by 1 (at least for > 1). Thus, 
this instability is not independent of the one analyzed in 
[25]. Since in our numerical work we focus on the Higgs 
sector, we will also focus our analytical considerations 
mostly on the instability mode discussed above (6), al- 
though we will estimate the energy gain and loss by this 
process at the end of this section. 

We will now show that in the presence of a background 
thermal bath of photons, the Z string is stabilized. For 
our analysis to hold, it is important that all of the scalar 
fields be out of thermal equilibrium, that the photon is 
in thermal equilibrium, and that the other gauge fields 
be out of equilibrium. These conditions are naturally 
satisfied in the electroweak theory at temperatures below 
the electroweak symmetry breaking scale and below the 
mass of the Higgs particle, but above the temperature of 
recombination. In this (large) temperature interval it is 
justified to average over the light degrees of freedom in 
order to study the dynamics of the order parameter. 

Thus, our procedure will be to consider the Lagrangian 
of the electroweak theory in the presence of a thermal 
bath of photons. We take the thermal average of this La- 
grangian and extract the terms which act as a correction 
to the potential for the dynamics of the order parameter. 
The thermal averaging consists of setting averages of 
to zero, and making the replacement 

A^A^ ^ -aT^ , (9) 

where a is another positive constant of order 1 (the es- 
timate of [4] gives a value somewhat larger than 1, and 
the larger the number of degrees of freedom in thermal 
equilibrium in the plasma, the larger the value of a will 
be). 

Thus, the starting point is the Lagrangian (2) . In this 
Lagrangian. we set the charged gauge fields to zero, and 
invert the transformation (3) of the neutral fields in order 
to express the Lagrangian in terms of Af^ and Z^. If 
the charged scalar fields are excited as in (6), then the 
transformation (3) is effected. However, the efiiect will be 
a correction of order in a quantity which is already of 
order and can hence be neglected. Thus, we use 



= cos{0^)Z + sm{e,,)A 
Y = cos{0u,)A-sm{e,^)Z . (10) 

Inserting these transformations into (2) allows us to ex- 
tract the extra contribution to the potential energy den- 
sity of the scalar fields which stems from the presence of 
the bath of photons. This contribution comes from the 
part of the covariant derivative term which is quadratic 
in A. Using for ^ the Z string configuration, it is easy to 
verify that the effective potential becomes 

Veff =Vo + l< AiA' > (25sin(e^))'($+)t$+ 

= Vo + ia^2(2gsin(^^^))'($+)t$+ , (n) 

where Vq is the bare potential (the last term on the right 
hand side of (2)). 

From Eq. (11) we see that the interaction with the pho- 
ton plasma induces a positive contribution to the mass of 
the charged scalar doublet. For temperatures in excess of 
a new critical temperature T^, this positive contribution 
will be larger than the negative contribution to the mass 
(when expanded about $ = 0) from the bare potential 
Vq. For temperatures smaller than T^, the total mass 
term is negative. From (11) and (2) we can immediately 
read off the value of T^: 

/2A\^^^ 1 
"^"^ ^ \ " / fi'sin(e^) ■ ^^■^^ 

We thus expect the core of the embedded defect to be 
symmetric (no charged scalar fields excited) for T > T^, 
whereas we expect a core phase transition [8] (charged 
scalar fields non-vanishing) for T < Td- However, if we 
consider the time evolution of an embedded string which 
is initially set up with symmetric core (as we will do in 
the simulations described below), then the temperature 
Td at which the symmetric vortex undergoes a core phase 
transition is expected to be lower, since gradient energy 
is required in order to produce an asymmetric core, not 
just the potential energy which enters the above consid- 
erations. 

Let us return to the issue of W-condcnsatc formation 
and briefly discuss the stabilization mechanism from this 
point of view. The energy gain per unit string length 
(computed by integrating from the center of the string 
to radius p) at zero temperature obtained from the for- 
mation of a condensate in which the amplitude of the 
W-field at the core radius is denoted W is of the order 
[26] 

i?(p)gai„ ^ g^ri^W^p" (13) 

(this excludes the gradient energy required to form a core 
condensate). At finite temperature, there is an energy 
loss due to the interaction of the photon with the W- 
fields. This energy loss was used in [28] to show that Z 
strings in strong magnetic fields can be stabilized, and the 



3 



corresponding energy per unit string length (integrated 
from the center to radius p) is given by 

E{p\ss ~ e^aT^W^p^ . (14) 

As is obvious from comparing (13) and (14), at suffi- 
ciently high temperatures the electrowcak string will be 
stable towards the formation of a W-condensate. The 
critical temperature is of the order a~^/'^ri, as in (12). 
3. Core Phase Transition 

We have simulated the evolution of embedded defects 
in the presence of a finite temperature charged plasma 
using a numerical code based on the one employed in [1] . 
Although gauge fields arc not included in this case, the 
stability of the string configuration can be established 
since it essentially comes from the modification of the ef- 
fective potential for the scalar fields and the gauge field 
should play an important role only when one string in- 
teracts with another one. 

First we set up the initial configuration as a two- 
dimensional slice of an infinitely long straight global 
string which is formed by the neutral components of the 
scalar fields and whose center resides in the midpoint of 
the lattice. Although such a highly symmetric configu- 
ration might be too ideal, it would be appropriate to see 
whether the core phase transition or the complete decay 
of the string occurs since the conservation of the winding 
number and the scalar field structure at the string core 
can be easily checked in contrast to full three-dimensional 
simulations. We then add thermal energy to the config- 
uration in the form of kinetic energy, that is, the time 
derivative of the scalar fields. Its amplitude is 0.1 x T^, 
and the allocation to the four components of the scalar 
field is chosen at random. 

Then, the four scalar fields (j) are evolved numerically 
on a two-dimensional lattice by means of the equations 
of motion derived for the scalar field with the effective 
potential (11). Neumann boundary conditions are em- 
ployed. During each simulation, the backgroimd temper- 
ature is constant and the cosmic expansion is not taken 
into account. However, in order to reduce the fluctu- 
ations of the fields so that it is easier to see whether 
the string configuration is preserved or not, an artificial 
damping term is introduced. Thus, the evolution equa- 
tions can be written as 

— -Vc^ + D- = -—, (15) 

and the numerical value of D is set to be 0.1 in the simula- 
tion. The inclusion of this damping term and the choice 

of the numerical value of D do not affect the essential 
results concerning the stability of the string. All dimen- 
sional quantities are rescaled by appropriate powers of 
the symmetry breaking scale, rj and are made to be di- 
mensionless. Although most simulations arc performed 
in a box size of 1000^, we have checked that the basic 
results are insensitive to the box size by executing 3000^ 



and 6000^ grid point simulations. The spatial resolu- 
tion is Ax = 0.5??~^, and the time steps are chosen as 
At = ^Aa;. The values of the numerical parameters are 
chosen as A = 2.5 x 10~^ and a{2gsuY{9w)Y = 2. We 
have calculated various patterns of the temperature and 
some of the results are depicted in the figures. 

Figures 1 & 2 show the resulting scalar fields as a func- 
tion of time and averaged over various sizes of square 
boxes centered at the midpoint of the simulation box, 
that is, the string core of the initial configuration, for a 
high temperature of T = rj. Figure 1 depicts the am- 
plitude of the neutral scalar field averaged over boxes 
of 10^ and 50^ grid points as well as averaged over the 
entire volume (1000^ grid points). Figure 2 shows the 
average of the charged scalar field over the entire volume 
(the average values over the smaller volumes are almost 
identical). 

The width of the field configuration of a Niclsen-Olesen 
string is about A~^/^. Thus, the averaging in the small- 
est of our three volumes (the 10^ grid point volume) is 
probing mainly the core region of the embedded string, 
and thus the amplitude of the neutral field component is 
small compared to unity, which is equal to the symmetry 
breaking scale, ry, in our normalization scheme, whereas it 
is already of order unity on the scale of the 50^ grid point 
volume. The fact that the neutral scalar field averaged 
over the smallest box (enclosing the initial electroweak 
string core region) remains small indicates that the elec- 
troweak string does not decay. The fact that the charged 
scalar field averaged over the core region remains vanish- 
ingly small indicates that no core phase transition takes 
place: the defect is a symmetric embedded defect. 

Figures 3 fc 4 show the corresponding curves in the 
case of a much lower temperature T = IQ^^rj. Note that 
this temperature is lower than the critical temperature 
Trf when the effective barrier stabilizing = in the ef- 
fective potential VcS disappears. The field configuration 
begins in the same state as in the high temperature case. 
However, after a short time; a core phase transition sets 
in during which the charged scalar field components take 
on a non- vanishing value which is of the order of rj. The 
fact that the neutral scalar field remains small when av- 
eraged over the smallest of the three boxes demonstrates 
that the embedded string remains stable. 

In fact, the absolute value of the neutral scalar field 
components is observed to decrease slightly during the 
core phase transition and, in contrast, the amplitude of 
the charged scalar fields averaged over the whole sim- 
ulation box is increasing. This represents the gradual 
increase in the core size, an increase which will stop once 
the core size R is of the order T~^. This can be seen as 
follows: If the core phase transition occurs via the ex- 
citation of a single of the two charged scalar fields (no 
winding number in the charged scalar field sector gen- 
erated), then the energy is lowered by eliminating the 
angular gradient energy of the neutral scalar fields. The 
energy per unit length thus gained will be of the order 
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FIG. 1. The amplitude of the neutral scalar field, 



averaged over volumes of 10 (dashed line) and 



50 (dotted line) grid points, and over the entire volume(solid 
line), as a function of time, in the high temperature simulation 
with T = 7?. 




1000 2000 3000 
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FIG. 3. The amplitude of the neutral scalar field averaged 
over volumes of 10^ (dashed line) and 50^ (dotted line) grid 
points, and over the entire volume(solid line), as a function 
of time, in the low temperature simulation with T = 10~^?7. 
The core phase transition occurs at a time t ~ 800 and is 
marked by a sharp docroaso in the average value of field in 
the string core and surrounding regions. 
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FIG. 2. The amplitude of the charged scalar field, 
\/ <jyi + 4>2^ averaged over the entire volume, as a function of 
time, in the high temperature simulation with T = rj. 
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time 

FIG. 4. The amplitude of the charged scalar field averaged 

over volumes of 10'^ (dashed line) and 50^ (dotted line) grid 
points, and over the entire volume(solid line), as a function 
of time, in the low temperature simulation with T = 10~^»7. 
At the time of the core phase transition t ~ 800, the charged 
scalar field takes on an average value over the string core 
region and the neighboring volume of 50^ grid points which 
is of the order r}. 
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R 



(16) 



where w is the core width before the core phase transi- 
tion. However, there is an energy cost associated with 
the generation of a nonvanishing value for the charged 
field, and this energy cost (per unit string length) is pro- 
portional to 



(17) 



By balancing the energy gain and the energy loss, one 
obtains an optimal core width which is of the order of 
R ^ T^^. Based on this consideration, we predict that 
the increase of the average value of the charged field over 
the entire box will eventually come to a halt, on a time 
scale which is proportional to T~^. We have confirmed 
these predictions concerning the final value of R and the 
time scale in our numerical work. 

The time at which the core phase transition takes place 
depends sensitively on the value of the damping term 
added to the equation of motion; the smaller the damp- 
ing effect is, the earlier the transition occurs. However, 
such a difference is not cosmologically significant when 
we consider the stability of the electroweak string since 
the time it takes for the core phase transition to occur 
is much smaller than the Hubble time both in either the 
cases of T = lQ~^ri and T = rj. The asymptotic am- 
plitudes of the fields averaged over the larger two of the 
three volumes do not depend on the presence or absence 
of the damping term, although the amplitude of the fields 
in the smallest volume does. This is due to the fact that 
in the absence of damping, the thermal fiuctuations will 
cause the string core to move by more than a core radius. 

Note that in both the high and the low temperature 
simulations, the winding number (in the neutral scalar 
field sector) is conserved not only for the entire simulation 
box but also when it is calculated around the string core, 
the 10^ grid points square region. This is a further test 
of the stability of the embedded string. 

We have also investigated numerically at which tem- 
perature Td the core phase transition takes place. We 
expect [1] the phase transition to happen at the temper- 
ature Td when the potential barrier at (/) = due to the 
plasma terms disappears. This occurs when the positive 
quadratic contribution to V^ofr due to the plasma equals 
the negative quadratic contribution due to the bare po- 
tential (see (12)). For the values of A and a{2g sm{9m))^ 
chosen the value of T4 is O.I77. However, an initially sym- 
metric core can only undergo a core phase transition if 
there is sufficient energy released from potential energy 
to create the required gradient energy. Our numerical 
simulations show that this happens at a temperature Td 
which lies in the interval between 0.04 and 0.05 in units 
of r] (see Figures 5 & 6). Note that the numerical value 
of Td for diff'erent model parameters such as A or a can 
be obtained by simple scaling using (12). 

Moreover, we have investigated the stability of the 
embedded electroweak string at very low temperatures, 
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FIG. 5. The amplitude of the neutral scalar field aver- 
aged over the entire volume, as a function of time, for 
four different temperatures near the critical temperature Td, 
T = O.02?7(long dashed line), T = 0.03r; (short dashed line), 
T = 0.047?(dotted Une) and T = 0.057?(solid line). 
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FIG. 6. The amplitude of the charged scalar field av- 
eraged over the entire volume, as a function of time, for 
four different temperatures near the critical temperature Td, 
T = O.0277(long dashed line), T = O.OSt? (short dashed line), 
T = 0.04?7(dotted Une) and T = 0.05??(soUd line). 
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T = 10~^rj. By tracking the winding number in the neu- 
tral field sector, as well as by tracking the ratio of the 
neutral to the charged scalar fields in the string core re- 
gion, we find no evidence for a decay of the string. The 
time evolution of the fields is almost as shown in Figs. 3 
and 4, except for the fact that the core phase transition 
sets in later since the thermal fluctuations are weaker. 
Obviously, our approximate analysis is only valid at tem- 
peratures for which the photon is in equilibrium, and thus 
breaks down before the cosmic temperature reaches that 
of last scattering. 

4. Conclusions and Discussion 

We have studied the stability of the embedded Z-string 
of the standard electroweak theory in the presence of an 
electro-magnetic plasma in an approximate treatment in 
which we only follow the dynamics of the Higgs fields 
and treat the gauge fields as either vanishing (the W 
and Z fields) or else (the photon) as being in a thermal 
bath. We find that the electroweak string is stabilized at 
all temperatures above which the photon is in thermal 
equilibrium. In the temperature range < T < 
the embedded string is symmetric in the sense that the 
charged scalar fields are not excited in the core, for lower 
temperatures the charged scalar fields arc non- vanishing 
in the string core, but the string itself is preserved in the 
sense that the winding number in the neutral scalar field 
sector is conserved. 

The basic mechanism which stabilizes the electroweak 
Z-string is the plasma mass for the charged Higgs doublet 
induced via the interactions with the plasma. This lifts 
the vacuum manifold in the direction of the charged Higgs 
doublet, leaving an effective vacuum manifold = iS^ 
which admits stable cosmic string solutions * 

Our numerical work is based on simulations in which 
only the scalar fields are evolved. From a gauge theory 
perspective this is an inconsistent approach. However, 
we feel that since the stability of the embedded defects 
is determined by the scalar field effective potential, ne- 
glecting the gauge fields should not adversely affect our 
results. However, it would be interesting to perform full 
local field theory numerical simulations to verify our con- 
clusions. 

The plasma stabilization mechanism discussed in this 
work obviously extends to a wide class of gauge theories 
which might be relevant in the early Universe. A new 
class of stabilized topological defects will thus create a 
new and interesting avenue to explore the interface of 
particle physics and cosmology. 
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